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1 Introduction 



Einstein's equations have, for much of the history of general relativity, been 
explored very fruitfully in terms of their concise and elegant statements char- 
acterizing the geometry of four dimensional pseudo-riemannian geometries. 
Such geometries depict possible physical spacetimes containing only "the 
gravitational field itself." The variety and properties of these "empty" space- 
times is truly astonishing. Quasi-local geometric entities such as trapped sur- 
faces and event horizons have become familiar. It is now firmly established 
that large-scale topological and geometrical features of spacetime are, indeed, 
subjects of physical inquiry. The nature and distribution of "matter" at stel- 
lar scales and upward has also brought particle physics and hydrodynamics 
to the fore. 

During these years, however, a steady development of the "space-plus- 
time" or 3 + 1 view of spacetime geometry has also occured. Here one views 
general relativity as "geometrodynamics" in the parlance of John Wheeler |I| . 
The emphasis, in the canonical or Hamiltonian explication of geometrody- 
namics given by Arnowitt, Deser, and Misner ("ADM") and by Dirac || 
U, is on the evolving intrinsic and extrinsic geometry of spacelike hyper- 
surfaces which determine, by knowledge of the appropriate initial data and 
by classical causality, the spacetime "ahead" (and "behind"), if spacetime is 
globally hyperbolic, an assumption we adopt throughout. 

Underlying and preceding geometrodynamics and Hamiltonian methods, 
however, was the basic realization that four of the ten Einstein vacuum equa- 
tions are nonlinear constraints on the initial Cauchy data, which play such 
a decisive role in defining the later canonical formalism ||. The Cauchy 
problem, constraints plus evolution, was shown to be well-posed in the mod- 
ern sense of nonlinear partial differential equations || [5], §, ||. This train 
of progress was marked by early work of Darmois and Lichnerowicz M , and 
brought to definitive development by one of us || |7|. 

At this writing, with accurate three-dimensional simulations using the full 
Einstein equations, with and without the presence of stress-energy sources, 
becoming essential for realistic studies of gravity waves, high energy astro- 
physics, and early cosmology, studies of Einstein's equations of evolution in 
3 + 1 form have blossomed [11-28] (an incomplete sample — see also [|10|1 ). 
Hyperbolic forms, especially first-order symmetrizable forms possessing only 
physically causal directions of propagation, have undergone very significant 
development in the past few years. 
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Specifically, in this chapter we describe in detail several explicitly causal 
hyperbolic formulations of Einstein's dynamical 3 + 1 equations by following 
a path that can be viewed as lighted by the "slicing function," a — the 
quantity that relates the lapse N to the determinant of the spatial metric g 
through N = g l ^ 2 a. This representation of the lapse function was presented 
in [Tj|. The slicing function allows us to demonstrate explicitly that no 
foliation of spacetime by spatial time-slices can be an obstacle to the causal 
hyperbolic forms of the dynamical Einstein equations. The slicing function 
plays an essential role (1) in a more precise form of the canonical action 
principle and canonical dynamics for gravity, (2) leads to a recasting of the 
Bianchi identities VpG^a = as a well-posed system for the evolution of the 
gravitational constraints in vacuum and also (3) of VpT^a = as a well-posed 
system for evolution of the energy and momentum components of the stress 
tensor in the presence of matter, (4) in an explicit display of four hyperbolic 
formulations of Einstein's equations with only physical characteristics, and 
(5) even in providing guidance to a new elliptic "conformal thin sandwich" 
form of the initial value constraints. 

We recall that the proof of the existence of a causal evolution in local 
Sobolev spaces of g and its extrinsic curvature (second fundamental tensor) 
K into an Einsteinian spacetime does not result directly from the equations 
giving the time derivatives of g and K in terms of space derivatives of these 
quantities in a straightforward 3 + 1 decomposition of the Ricci tensor of 
the spacetime metric, which contains also the lapse and shift characterizing 
the time lines. These equations do not appear as a hyperbolic system for 
arbitrary lapse and shift, in spite of the fact that their characteristics are 



only the light cone and the time axis [25 



We now turn to notational matters, conventions, and to the 3 + 1 decom- 
position of the Riemann and Ricci tensors. We assume here and throughout 
the sequel that the spacetime V = M x K is endowed with a metric g of 
signature (— , +, +, +) and that the time slices are spacelike, that is, have 
signature (+,+,+). These assumptions are not restrictive for globally hy- 
perbolic (pseudo-riemannian) spacetimes. 

We choose on V a moving coframe such that the dual vector frame has 
a time axis orthogonal to the slices M t while the space axes are tangent to 
them. Specifically, we set 

6° = dt, 

e i = + (1) 
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with t 6 M and x l , % = 1, 2, 3 local coordinates on M. The Pfaff or convective 
derivatives d a with respect to 8 a are 

In this coframe, the metric g reads 

ds 2 = g a p6 a 6P = -N 2 (6 ) 2 + g ij 9 i 9 j . (3) 

The t-dependent scalar N and space vector (3 are called the lapse function 
and shift vector of the slicing. These quantities were explicitly identified 
in [0 and play prominent roles in all subsequent 3 + 1 formulations. Any 
spacetime tensor decomposes into sets of time dependent space tensors by 
projections on the tangent space or the normal to M t . 

We define for any t-dependent space tensor T another such tensor of the 
same type, d T, by setting 

do = | - £p , (4) 

where £p is the Lie derivative on M t with respect to (3. 

Notice that in our foliation-adapted basis (|I|) and (0), that if g denotes 
the spacetime metric and g the space metric, then we have (g 0i = g 0t = in 
our frames): 

9ij = 9ij ; g 13 = g 13 ■ (5) 

(Greek indices range {0, 1,2,3} while latin ones are purely spatial.) Hence, 
no overbars will be used to denote components of the spatial metric. On the 
other hand, for the determinants, we have (—detg) = N 2 (detg). Therefore 
we shall use overbars on spatial metric determinants; for example g 1 ! 2 = 
(detg) 1 / 2 . 

Likewise, to distinguish the purely spatial components of the spacetime 
Ricci tensor (say), we shall write Rij(g), while for the space Ricci tensor we 
shall write Rij(g). In general, of course, Rij(g) ^ Rij(g). The Levi-Civita 
connection of g is denoted by V and that of g by V. 

With the convention 

\7 a ap = d a ap - a p Y/3a , (6) 
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and the definitions 

7°/3 7 = + 9 a& C 6 8(^)6 - Tp a Pl ' ( 7 ) 

do a = -\c a Pl e^Ae\ (8) 

we have for the connection coefficients (r denotes an ordinary Christoffel 
symbol) 

Yj k = r'Ag) = r'Ag) (9) 

Yo k = -NK\ , 7 l ,o = -NK\ + dj? , 7 ° - = -iV" 1 ^ (10) 

7*00 = N&N , 7°oi = 7°<o = 9ilogiV, 7 °oo = <9 log iV . (11) 
Observe that if a is a space scalar of weight —1, we have 

V,a = + ar fc fci ((?) = d z a + adi log^ 1 / 2 , (12) 

£ pa = ftVia + aVifr . (13) 

The Riemann tensor is fixed by 

(V^ - V^V a ) ]T = V 5 W SaP (14) 

while the Ricci tensor is Rsp = -R 7 «5 7/ 3- 

The 3 + 1 decompositions of the Riemann and Ricci tensors are 

Rijki(g) = Rijki(g) + 2K i[k K i]j , (15) 

R 0ijk (g) = 2NV {j K k]l , (16) 

RoM = N(d K lJ +NK lk K k J + V t d J N) . (17) 

One can then obtain for the Ricci tensor 

R tJ (g) = R ij (g)-N- 1 d K ij + KK ij -2K ik K k j -N- 1 V i d j N , (18) 

R 03 (g) = N{d j K-V h K h ] ) , (19) 

R 00 (g) = N{d K-NK ij K i ^A- g N) , (20) 

where K = K% and A 5 = g^VjVj. Finally, we note 

G\ = \(K l3 K^-K*-R{g)) . (21) 
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2 Every Time Slicing Is "Harmonic" 

The standard statement of the harmonic time-slicing condition is, that on a 
t = const, time slice, Bq[(— g) l ^ 2 g m ] = 0. (This is equivalent, in a coordinate 
basis, to dfj[{— gY^g^] = 0.) Friedrich observed in |24| that the right hand 
side of these equations could be a given function of (t,x l )Q. (See also J0J.) 
Therefore, the standard harmonic condition expressed in 3 + 1 form, d$N + 
N 2 K = 0, can be written as a generalized "harmonic" condition 

B N + N 2 K = Nf , (22) 

where f(t,x) is a known function. Specifically, introduce a(x, t) such that 
<9ologa = /, then (p2|) becomes 

d N + N 2 K = Nd log a , (23) 

from which the identity 

Bologg 1 / 2 = -NK , (24) 

allows us to see that 

N = g 1/2 a . (25) 

We shall call a(x,t) the "slicing function;" it is a freely given scalar density 
of weight —1. 

It is clear that any iV > on a given time slice t = to can be written 

1/2 

in the form N to = g t £ a(to,x) for some a > provided that gij(t ,x) is 
a proper riemannian metric. Introducing "harmonic" time-slicing is thus a 
simple matter. It is not, however, known at present how to construct a 
specific long-time foliation from general rules telling how to specify a(t,x). 
However, many foliations can be constructed in a "step-by step" fashion 
(numerical time steps) provided certain obvious conditions are met. For 
example, an elliptic condition on N can determine a(t) on a sequence of 
time slices if the condition does not couple to variables that disturb the 
characteristic directions of the hyperbolic equations. (The same is true for 
the shift vector (3 l .) Alternatively, we can try educated guesses for a(t,x). 

As it stands, (^3|) is clearly a speed zero (with respect to do) hyperbolic 
equation. However, this equation and Einstein's equations lead to a second 

1 Just as in electrodynamics, V^A^ = — * V^A^' — £(t,x) ^ is perfectly acceptable 
as a "Lorentz gauge" if £ is known. 
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order equation in space and time that propagates N along the light cone. This 
result brings into sharp relief the congruence of (p3|) with the propagation on 
the light cone of other variables. 



The trace of Rij(g) gives an equation for 8qK |L9 



B K = -A- g N + [R(g) +K 2 - R k k (g)} N , (26) 



where in ( |26|) denotes the Laplacian <7 y VjV 3 -. Taking the time derivative 
of ( p3|) and eliminating d$K with (^) shows that N obeys the non-linear 
wave equation 

U g N + R\{g)N - R(g)N - Nd log a + (d 2 log a)^ 1 = , (27) 

where we wrote our wave operator or "d'Alembertian" as □ = — (N~ 1 do) 2 + 
Ag. The characteristic cone of □ is clearly the physical light cone (c = 1). 
The equation (27) per se will not be used explicitly in the sequel. 

Substitutions of the form N\ = g x l 2 a\ (A > 0) have also been consid- 
ered p7| . However, after working out the wave equation analogous to (|27|) 



that A^a obeys, one finds that the local proper propagation speed of N\ is 
VX This behavior may or may not spoil the propagation of system variables 
other than N, but if A ^ 1 and the system is hyperbolic, one will always 
find that the characteristic directions of the system will not all be physical 
ones. That is, in vacuum gravity, there will be some variables that propagate 
neither on the light cone (speed = 1) nor along the axis parallel to do (speed 
= 0) which is orthogonal to t = const. The variables not propagating in 
physical directions are gauge variables and one will not have physical crite- 
ria for their boundary values on characteristic surfaces. On the other hand, 
with A = 1, one has fulfilled a necessary condition that physical and gauge 
variables propagate together in the same directions. 

In the following sections, whenever we consider hyperbolic systems, we 
will focus on first-order symmetric (or symmetrizable) hyperbolic ("FOSH") 
equations possessing only physical characteristic directions. We understand 
"FOSH" in this restricted physical sense only in this paper, and likewise for 
other uses of the term "hyperbolic." 

3 Canonical Action and Equations of Motion 



Choice of the slicing function a in ( pop is arbitrary (a > 0). That a is 
freely chosen while N must satisfy an equation of motion (|23|) suggests that 
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it should be regarded as the undetermined^ multiplier in the canonical action 



principle of Arnowitt, Deser, and Misner ("ADM") |2j. Here we follow |2S 



To draw some lessons for the canonical formalism, let us first express the 
3 + 1 evolution equations in their standard geometrical form (see: with zero 
shift [§], arbitrary lapse and shift @, spacetime perspective pOj): 

9ij = —2NKij , (28) 
Kij = N (-RM + Rij(g) + NKij - K ik K k - N^VftN) , (29) 
where ( ) = Bq( )■ 

A brief look at ( p9| ) shows that forming the combination TZij = Rij(g) — 
gijR k k{g) leads to an equation of motion for the ADM canonical momentum 

7r« = gW (Kg ij - K ij ) (30) 

that contains no constraints. (In this section we choose units in which 
167rG = c = 1.) Indeed, using (p8| ) and (]29|) , we obtain the identity 

tfi = Ng 1/2 (R{g)g lj - R ij (g)) - Ng' 1 ' 2 {2ix ik Ti j k - Tm ij ) 

+g 1/2 (V l V j N - g ij V k V k N) + Ng 1/2 [K ij ] . (31) 

From the identity gij = —2NKij, we have 

g .. = N g-W ( 27rfi - ngij ) . (32) 

We now come to a crucial observation. Were the canonical equation for 
jr^ to be dictated by vanishing of the spatial part of the Einstein tensor, 
G l i{g) = 0, as it is in the conventional ADM analysis 0, then the identity 

GM + g l0 G%{g) = R^g) - 9ij R k k (g) = 7^- (33) 

shows that a Hamiltonian constraint term ~ g 1 ^ 2 G°Q remains in the tt^ equa- 
tion (|3T|). This would mean that the validity of the tc^ equation would be 
restricted to the subspace on which the Hamiltonian constraint is satisfied 
(i.e., vanishes). 

2 The multipliers associated with gauge freedom or, as here, with spacetime coordinate 
freedom, can be freely chosen because they are not determined by physical conditions. 
Hence, they are not true "Lagrange multipliers." 
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Though the ADM derivation of the ir^ equation, found by varying in 
their canonical action (f3 l is the shift vector) 

S [g, vr; N,(3) = J d A x (tt% - NH) , (34) 

with N(t,x) , /3 l (t,x) and 7r u held fixed, is of course perfectly correct, an- 
other point of view is possible. [We are ignoring boundary terms, a subject 
not of interest here, and we note that the momentum constraint term — (3 l 7ii 
(Hi = g 1//2 Cj, Ci = 2NR°i) is contained in n^gij (( ) = 8q) upon integration 
by parts.] (The slicing density a has also been used prominently in the action 
by Teitelboim |TT], who simply set a = 1 (N = g 1 ^ 2 ), and by Ashtekar 
for other purposes.) 

We have explained that a can be regarded as a free undetermined multi- 
plier while N is a dynamical variable (a conclusion also reached by Ashtekar 
for other reasons p2| , [331 ) that determines the proper time N5t between slices 
t — t' and t = t' + 5t. N is determined from a(t, x) and g 1 ^ 2 found by solving 
the initial value constraint equations. (See the treatment of the constraints 



in the final section of this article and in |sy, |34|, |35|, [|.) Motivated by this 
viewpoint, we alter the undetermined multiplier N in the ADM action prin- 
ciple to a, where the Hamiltonian density H is (with H. = 2g 1 / 2 G°o(g) being 
the ADM Hamiltonian density of weight +1) 

U = g 1/2 H = - \* 2 - gR(g) , (35) 

a scalar density of weight +2 and a rational function of the metric. The 
action becomes 

S [g, it; a, (3) = J d 4 x (n^cjij — afi^j . (36) 

The modified action principle for the canonical equations that we propose 
in fl36l) is to vary ^ and g^j, with a(t, x) and /3 l (t,x) as fixed undetermined 
multipliers. From 

8H = (2^-^vr) Sir" + (27r ik ni k - + gRV(g)-ggVR(g))5 gij 

-^(W^y-^VfcV*^) , (37) 

we obtain the canonical equations 

9v = a J~Tj = a ( 2n ij ~ n 9ij) = -2NK i:j , (38) 
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tt'J = -a^- = -ag (R ij (g) - R(g)g ij ) - a (2n ik ^ k - 7nr ij ) 

+g(V l V j a-g ij V k V k a) . (39) 
Equation (|39|) for n 1 ^ is the identity (|3~TD with IZ 1 ^ = 0, which is equivalent to 



Rij(g) = 0. Thus, ( p9|) is a "strong" equation unlike its ADM counterpart, 
which requires in addition the imposition of a constraint: TC — 0. 

In the present formulation, the canonical equations of motion hold every- 
where on phase space with any parameter time t, a necessary condition for 
the issue of "constraint evolution" even to be discussed in the Hamiltonian 
framework. (See below in Sect. [![) 

If we define the "smeared" Hamiltonian as the integral of the Hamiltonian 
density, 

H a = d 3 x'a{t,x')H , (40) 



the equation of motion for a general functional F[g, 7r; t, x) anywhere on the 
phase space is 

F[g,Tr;t,x) = -{n a ,F^+d F , (41) 

where ( ) denotes our total time derivative and do is a "partial" derivative of 
the form d t — £p acting only on explicit spacetime dependence. The Poisson 
bracket is 

{F G\=!d* X (^ ^ ^ 5 -^—\ (42) 

J \5g ij (t 1 x)5ir i i(t,x) 5 9ij (t,x) 5^(t,x) J ' v ' 

and one sees that time evolution is generated by the Hamiltonian vector field 

X na = J d 3 x |a(2 7 r ij - ^)^- - [ag(R ij (g) - R(g)g ij ) 

+ a (2n lk nh-rrn^)~g(V l V j a-g l ^ k V k a)]^--^ . (43) 

Because it does not contain any explicit constraint dependence, ([43D is a valid 
time evolution operator on the entire phase space. It is clear that the (g, 7r) 
equations come from ( f4~3D applied to the canonical variables. The harmonic 
time slicing equation ( |23| ) results from application of ( f4"3"D to N, and the wave 
equation for iV comes from a repeated application of (f43D to (|23|). 
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Evolution equations for the "constraints" are computed to be 



d H = -{H a ,H} =agg^d l H J + 2gg i m t V J a , (44) 
BoHj = - {W a , Wj} = aSjW + 2W0 i a , (45) 

where Vj-a = o^a + ag~ l ^ 2 dig 1 / 2 . These are well-posed evolution equations 
for the constraints, and they are equivalent to the twice-contracted Bianchi 
identities when TZy = or Rij = (see below). 

These results shed new light on the Dirac "algebra" of constraints 



It is well known that the Dirac algebra is not the spacetime diffeomorphism 
algebra. This can be seen from the fact that while the action (RBI) is invariant 



under transformations generated by TCj and TC, |37| the equations of motion 
that follow from this action are Rij{g) = even when Tij and 7i do not 
vanish. These equations of motion are preserved by spatial diffeomorphisms 
and time translations along their flow in phase space, whereas a general 
spacetime diffeomorphism applied to Rij(g) = would mix in the constraints. 

A second important view of the Dirac algebra results from the direct 
and beautiful dynamical meaning of its once-smeared form. Equations (53) 



and ( f45|) express consistency of the constraints as a well posed initial-value 
problem. If the constraint functions vanish in some region on an intial time 
slice, they continue to do so under evolution by the Hamiltonian vector field 
into the domain of dependence of that initial region. This mechanism follows 
from the dual role of 7i as a constraint and as part of the generator of time 
translations of functionals of the canonical variables anywhere on the phase 
space. 

Let us take note that the Hamiltonian constraint per se does not express 



the dynamics of the theory; the equation of dynamics is (f|l]). In its "altered" 
role, the Hamiltonian constraint function simply vanishes as an initial value 
condition , from which g 1 ' 2 is determined as in the initial value problem. BO] 



Then N can be constructed from a. The Hamiltonian constraint, once solved, 
remains so according to the results embodied in (HI) and 



4 Contracted Bianchi Identities 

The results on canonical dynamics that follow on using a as an undetermined 
multiplier are also reflected in the manner in which the twice-contracted 
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Bianchi identities, 

VpG p a = , (46) 
can be written as a first-order symmetrizable hyperbolic system ]29|]. (In the 



absence of hyperbolic form, (fig ) is practically useless in providing physical 
equations of motion for the constraints when they are not satisfied.) Likewise, 
this system extends to matter (see below). (Frittelli obtained well-posedness 
for fl4"6j ) by other methods [0.) 

We recall that the equations of motion of the canonical momenta in vac- 
uum are 

K ij = R ij {g)-g ij R\(g) = U, (47) 
while the weight zero Hamiltonian constraint is 

C = 2G° (g) = K i3 K^ - K 2 - R(g) = , (48) 

and the weight zero one-form momentum constraint is 

d = 2NR°i(g) = 2V\K l3 - K 9ij ) = . (49) 

Recall the identity (|33| ): 

Gij(g) + gijG\{g) = R v (g) - g l3 R k k (g) = H tJ . (50) 

Combining (FT?]), (|49"D, and (P3|) with (EH) gives the twice-contracted 



Bianchi identities clS db FOSH system 

C - NV J Cj = 2 (Cj V j N + NKC - NK ij [K^]) , (51) 

C j - NVjC = 2 (cVjN + l -NKC 3 - V l (iV[7^])) . (52) 

Substituting Hi = g x l 2 C^ H = gC, and setting the equations of motion 
Hij = in (fy]) and ([52]) yields the evolution equations of the unsmeared 
constraints as in (fPj) and (|5|). 

Similar considerations show how to put the "matter conservation" equa- 
tions V/jT ^ = into well-posed form. This was also carried out by one of 
us (YCB) and Noutchegueme |$9| but the results obtained here are more im- 
mediately physical. Unlike P5| , we use the energy density e = — T° rather 
than p 00 (p a /3 = T a j3 — ij^T 7 ^) to obtain this result. It is clear that such a 
result is possible because 

H af3 = K~ l G aj3 (g) - T a/3 (53) 
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vanishes as Einstein's equation and in any case satisfies V p H p a = 0. We 
can treat H a p as we did G a p above, (k = 8ttG ; c = 1.) The result is 
nevertheless of interest as it presents the continuity and relativistic Euler 
equations of matter in a well-posed form. 

Straightforwardly expanding V/jT^o = and V^T^j = gives the conti- 
nuity and Euler equations (cf. |3H| , p. 89), with e = — T° and the matter 
current one-form ji = NT°i. The continuity equation is 

B e + NV'ji = N(K i:j T ij + Ke - 2j i( ?) , (54) 

where a, = V$ log N is the acceleration of observers at rest in a given time- 
slice. Likewise, we find for Euler's equation 

B ji + NVjT j i = Mix 1 ; - T'ia, - :a ; ) . (55) 



The divergence term on the left side of ( pq) spoils the well-posed FOSH 
form we seek. However, if we use the identity Gij(g) + gijG° (g) = Rij(g) — 
9ijR k k(g) and the Einstein equations K _1 G Q/ g — T a/3 = 0, or K~ l R a p = p a p, 
we obtain (e = -T° , j t = NT°i) 

T\ - 6\e = {(9, - S j iP k k ) = Sh . (56) 

Then ( |54"D and (55) obtain well-posed form (if is assumed known), 

B e + NV*ji = -2jiV*N + 2NKe + NK ij S i:j , (57) 

B oJt + NViS = -2eV t N + NK 3l - V'(JVSy) . (58) 

By combining ( |57D plus (|49"D, and ( |58"D plus fl50D, we obtain expressions 
of gravity constraint evolution in the presence of matter, 

B C T - N^Cf = 2 (CjV j N + NKC T - NK^^U^ - <%)) , (59) 



B Cj - NV C T 



C T V 3 N + ^NKCj - V 4 (NiK^TZij - <%)) 



(60) 



where C T = C + 2e and Cj = Cj — 2jj. This is just the form we would 
anticipate on the basis of Hamiltonian dynamics and the form ( |4"9"D and ( 50 ) 
of the vacuum constraints. Thus, for gravity plus a matter field, we obtain 
results analogous to (§3]) and ( |45"D for the total system. If there are no 
violations of constraints, then C T = , Cj = 0, while the dynamical gravity 
equation is K^TZy — = 0. 
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5 Wave Equation for Kij 

Einstein's equations, viewed mathematically as a system of second-order par- 
tial differential equations for the metric, do not form a hyperbolic system 
without modification and are not manifestly well-posed, though, of course, 
physical information does propagate at the speed of light. A well-posed hy- 
perbolic system admits unique solutions depending continuously on the initial 
data and seems to be required for robust, stable numerical integration and for 
full treatment by the methods of modern analysis, for example, exploitation 
of energy estimates. The well-known traditional approach achieves hyper- 
bolicity through special coordinate choices.^ The formulation described here 
permits coordinate gauge freedom. Because these exact nonlinear theories 
incorporate the constraints, they are natural starting points for developing 
gauge-invariant perturbation theory. 

Consider a globally hyperbolic manifold V = £ x M. with the metric as 
given in the introduction. To achieve hyperbolicity for the 3 + 1 equations, 
we proceed as follows. 

By taking a time derivative of Rij(g) and subtracting appropriate spatial 
covariant derivatives of the momentum constraints, one of us (YCB) and 
T. Ruggeri |T8| (see also |T9|] , where the shift is not set to zero) obtained an 



equation with a wave operator acting on the extrinsic curvature. In vacuum, 
one finds 

doRijig) - V.Rqj - VjR i = NU^ + J tj + Sy = , (61) 

where □ = — (iV -1 ^) 2 + VfcV fe , consists of terms at most first order 
in derivatives of K^, second order in derivatives of gij, and second order in 
derivatives of N, and 

Sij = -N^ViVjidoN + N 2 K) . (62) 

The term is second order in derivatives of and would spoil hyper- 
bolicity of the wave operator □ acting on K^. Hyperbolicity is achieved by 
setting N = g x l 2 a{t, x), or 

B N + N 2 K = g 1/2 B a(t, x) , (63) 



3 The classic second-order fully harmonic form was given in || [7| and discussed, for 
example, in ||. It will not be discussed in this article. A FOSH form based on these 
equations was given first by Fischer and Marsden in f23fl. 
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as discussed in Sect. 0. The resulting equation combined with ( B5J ) forms 
a quasi-diagonal hyperbolic system for the metric with principal opera- 
tor doO. This system can also be put in first order symmetric hyperbolic 



form [pll 11], by the introduction of sufficient auxiliary variables and by use 



of the equation for i?oo (thus incorporating the Hamiltonian constraint ). 
The Cauchy data for the system (in vacuum) [18, IS] are (1) (g,K) such 
that the constraints Rqh = 0, G°o = hold on the initial slice; (2) d^K^ 
such that Rij = on the intial slice; and (3) N > arbitrary on the initial 
slice. Note that the shift ft k (x,t) is arbitrary. Using the Bianchi identities, 



one can prove (18], [L9| that this system is fully equivalent to the Einstein 
equations. The point is that quasi- diagonal Leray hyperbolic systems 
have well posed Cauchy problems and therefore unique solutions for given 
initial data. Because every solution of the Einstein equations also satisfies 
the OKij equation in particular and provides initial data for it, uniqueness 
implies, conversely, that if the initial data for the OKy equation are Ein- 
steinian, all solutions of Einstein's equations, and only these, are captured. 
The restriction on the initial value of d^K^ prevents the higher derivative 
from introducing spurious unphysical solutions. 

All variables propagate either with characteristic speed zero or the speed 
of light. The only variables which propagate at the speed of light have the 
dimensions of curvature, and one sees that this is a theory of propagat- 
ing curvature. However, a FOSH system that propagates curvature is more 
transparent in the "Einstein-Bianchi" form (next section). 

In the above formulation, the shift and a are arbitrary. This and our other 
systems (except the Einstein- Christoffel system in Sect. 0) are manifestly 
spatially covariant and all time slicings (using a) are allowed. Spacetime 
covariance is therefore present, but not completely manifest. 

By taking another time derivative and adding an appropriate derivative 
of Rqo, one finds (in vacuum) p9[ 



dodoRij - d ViR 0j + doVjRoi + ViV 3 R 00 = d Q {NUK i3 ) + J i3 = , (64) 

where Ji 3 consists of terms at most third order in derivatives of gi 3 and 
second order in derivatives of Ki 3 . Together with d^g^, these form a system 
for (g, K) which is hyperbolic non-strict in the sense of Leray-Ohya. [O] Here, 
the lapse itself, as well as the shift, is arbitrary (N > 0). The Cauchy data 
of the previous form (in vacuum) must be supplemented by 8 d Kij such 
that doRij = on the initial slice. This guarantees that the system is fully 
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equivalent to Einstein's theory (except that its solutions are not in Sobolev 
spaces ||29|| ). This system does not have a first order symmetric hyperbolic 



formulation, but has been used very effectively in perturbation theory |42| 
and in other applications |^3|, |44 |. 



6 Einstein-Bianchi Hyperbolic System 

To obtain a first order symmetric hyperbolic system, one can use the Riemann 
tensor of the spacetime metric. It satisfies the Bianchi identities for the 
spacetime geometry 

Vo-R^A^t + V/3-R 7a AAt + V 7 -R Q , / 3A M = . (65) 

These identities imply by contraction and use of the symmetries of the Rie- 
mann tensor 

V QJ R> 7 + V 7 i^ + = . (66) 

If the Ricci tensor R a p satisfies the Einstein equations (k — c — 1) 

Ra/3 = Paf3 , (67) 

then the previous identities imply the equations 

V Q -R a ^/3 7 = V,ap 7M — V 7 p j a At . (68) 

The first equations with (a/3j) = (ijk) and the last one with /i = 
do not contain derivatives of the Riemann tensor transverse to M t . They 
are considered as "constraints" and will be identically satisfied (initially) in 
our method. They remain satisfied in an exact integration. All detail and 



rigor concerning this elegant system is given in |22|, [5T| , to which the reader 
is referred. It has 66 equations, just as do the Einstein-Ricci first order 
curvature equations. 



The system we are now developing |2T| is similar to an analogous system 
obtained by H. Friedrich [E^] that is based on the Weyl tensor. The Weyl 
tensor system is causal but with additional unphysical characteristics. 

We wish first to show that the remaining equations are, for n = 3 in 
the vacuum case, when g is given, a symmetric first order hyperbolic system 
for the double two-form R a p\^. For this purpose, following Bel Q we 



16 



introduce two pairs of "electric" and "magnetic" space tensors associated 
with a spacetime double two-form A, 



N 2 E l3 {g) = A 0i0j (69) 

DM = \e lhk e 3lm A hklm (70) 

NH i:i (g) = l -t lhk A h \ 3 (71) 

NB 3i {g) ee i,V fc e l/lfc (72) 



where e^fc is the volume form of g. It results from the symmetry of the 
Riemann tensor R with respect to its first and second pairs of indices (R is a 
"symmetric double two- form") that if A = R, then E and D are symmetric 
while Hij = Bji. A useful identity for a symmetric double two-form like R, 
with a tilde representing the spacetime double dual, is ( "Lanezos identity" ) 

Raf3\n + Raf3\n = C a \ g^ — C a ^ gp\ + Cp^ g a \ — Cpx g afl , (73) 

where C a p = R a p — (l/4)g a pR. It follows that when R a p = Xg a /3, then E = 
—D and H = B. In order to avoid introducing unphysical characteristics, 
and to be able to extend the treatment to the non-vacuum case, we do not 
use these properties in the evolution equations, but write them as a first order 
system for an arbitrary double two-form A, as follows: 

Vo^Ufcoj + VfcAohoj — V/jAofcoj = , (74) 
V A° i0j + V h A^ = V Pji - Vj p 0i , (75) 

and analogous equations with the pair (0 j) replaced by (Ira). One obtains a 
first order system for the unknowns E, H, D, and B by using the relations in- 
verse to the definitions above. The principal parts of these equations, all with 
one definite index fixed on E, H, D, and B, are identical to the correspond- 
ing Maxwell equations. The characteristic matrix of this "Maxwell" part of 
the system has determinant — iV 6 (£o£ ) (£ a £ a ) 2 - The system obtained has a 
principal matrix consisting of 6 identical 6 by 6 blocks around the diagonal, 
which are symmetrizable and hyperbolic. Hence, the system is symmetric 
hyperbolic, when g is a given metric such that g is properly riemannian and 
N > 0. 

To relate the Riemann tensor to the metric g we use the definition 

d o9ij = -2NK i3 - (76) 
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and we use the 3 + 1 identities given in the Introduction. Note that in this 
section all T 's are spatial. 

We next choose N = g l l 2 a(t, x). We generalize somewhat the ideas used 
by Friedrich (see ||25|| ) for the Weyl tensor to write a symmetric hyperbolic 



system for K and T, namely we obtain equations relating T and K, for a 
given double two-form A, and by considering the definition of K and the 
3 + 1 decomposition of the Riemann tensor, replacing in these identities the 
Riemann tensor by A. To deduce from this system a symmetric hyperbolic 
first order system, with the algebraic form of the harmonic gauge iV = g l l 2 a, 
one uses the fact that in this gauge one has 

T h ih = d i logN-d i ]oga. (77) 

We obtain 

d Q T\ 3 +NV h K l3 = NK ij9 hk {Y m mk + d k \oga) 

-IN&^^^ + dfiloga) (78) 
— N {t k {j h Bi)k + H k ^e k ,^ h ) , 

and 

BoKy + Nd h T\ 3 = N [T m ih T h jm - (T\ h + d t log a) (T k ]k + d 3 log a)} 

-N (d^ log a - Y k ij d k log a) (79) 
-N [D m + E m - D k k9ij - KKij] . 

The system obtained for K and V has a characteristic matrix composed of 6 
blocks around the diagonal, each block a 4 by 4 matrix that is symmetrizable 
hyperbolic, with characteristic polynomial — iV 4 (£ £°) {£, a £ a )- 

The whole system for A, K, T, g is symmetrizable hyperbolic, with char- 
acteristics the light cone and the normal to M t . It is somewhat involved to 
prove that a solution of the constructed system satisfies the Einstein equa- 
tions if the initial data satisfy the constraints, but we can argue as follows. 
We consider the vacuum case with initial data satisfying the Einstein con- 
straints. These initial data determine the initial values of T, and also, if f3 
and N are known at t — 0, the initial values of Aijhm, Ajhio, A i0 jh by using 
the decomposition formulas. (We set A equal to the Riemann tensor on the 
initial surface.) We use the Lanczos formula to determine Ajojo initially. We 
know that our symmetrizable hyperbolic system has one and only one solu- 
tion. Because a solution of Einstein's equations with iV = (7 1//2 a, proved to 
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exist in the section on □i^y, satisfies together with its Riemann tensor the 
present system and takes the same initial values, that solution coincides with 
the solution of the present system in their common domain of existence. 



7 Einstein-Christoffel System 

The first-order form of the wave equation for K^, the Einstein- Ricci sys- 



tem |r9| , |Tl|j , to which we have alluded, has 66 equations, the correct number 
for a curvature system as does the Einstein-Bianchi system of the previous 
section. It is symmetric hyperbolic and therefore well-posed. But it is natural 
to ask whether there is a simpler first-order system of fewer variables that is 
perhaps closer in form to (|28|) and (PS|) . Frittelli and Reula [^, |27| proposed 



such a system, but their system has unphysical characteristics and is not 
written fully in terms of geometric variables. Here we deduce a different sys- 
tem having only physical characteristics and expressed in geometric variables. 
We understand (private communication to JWY) that James Bardeen has 



likewise obtained a similar system improving that found in [|16|, [T7|. While 



our derivation |T(| can proceed systematically by direct construction of an 
energy norm and the characteristic speeds, a more heuristic derivation is in- 
dicated here from the structure of the wave equation for Ky. Not that in 
this section all T's are spatial. 

In the dynamical spacetime Ricci tensor Rij(g), (P9"D, one has a dynamical 
equation for the extrinsic curvature in terms of spatial derivatives of the 
spatial Christoffel symbols. When the lapse N is replaced by the slicing 
density through ag 1 ^ 2 , the differentiated Christoffel terms become 

dk^ k ij — djT k ik — diT k jk . (80) 

This may be read as the divergence of a linear combination of Christoffel 
symbols, which puts the Rij{g) equation in a form reminiscent of the structure 
of one of the first-order equations for a free wave, namely 

d u + d k v k = . (81) 

We would have a symmetric hyperbolic system if there were an analog of the 
other equation for a wave, 

d v k + d k u = 0. (82) 
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Some manipulation quickly leads to the conclusion that and the Christof- 
fel combination above are not paired in a symmetric hyperbolic system like 
u and Vk- 

Recall however, that the free wave equation (do) 2 u—d k dkU = is obtained 
by taking a time derivative of (|y]) and subtracting the divergence of (|32"|). 
In obtaining the wave equation for Kij (|6T|) we have taken a time derivative 
of Rij and subtracted a (suitably symmetrized) divergence of the momentum 
constraint Rq{. This motivates the speculation that in gravity the "other" 
equation should be related to the momentum constraint and its sole spatial 
derivative should be dkKij. 

Following this idea about the second equation leads one to consider an 
equation of the form 

dkiRjo + QkjRio = —d fkij — dk(NKij) + l.o.kij (83) 

where l.o.kij are lower order terms involving no derivatives of fkij or K^. 
One must choose fkij from a linear combination of spatial derivatives of the 
metric. Introduce 

Qkij = d k gij (84) 

and use the identity 

d {d kgij ) = -dkpNKij) (85) 

to find that 

fkij = Qphij ~ 9k{iQ TS {Q\rs\j) ~ Qj)rs) (86) 

produces the correct coordinate derivatives occuring in the momentum con- 
straints. The lower order terms are those terms necessary to complete ( j83|) 
into an identity and take the form 

l.O.kij = 2NK k (ig VS \G\rs\j) — Gj)rs) 

+2g k{i [Kj )m d m N -Kd 3) N 

+NK j)m g rs T m rs (g) + l -N (g, ya - 2Q lrslj) ) K rs \ , (87) 
where the spatial Christoffel symbols are constructed from Qijk, 

^kij(G) = (1/2) (Gjki + Qikj — Qkij) ■ (88) 

(It is clear that only one of the three-index symbols Qkij, Fkij, and fkij is 
needed, say fkij. The necessary algebra will not be reproduced here.) 
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One then easily verifies that by expressing flSS|) in terms of derivatives of 
the metric (assuming a metric compatible connection), we can manipulate it 
to take the form of a divergence of fkij plus lower order terms. The dynamical 
Ricci equation becomes 



where 



Rij = -N^doKij - d k f kij + l.o.ij , (89) 

l.o.ij = KK ij -2K ik K h j -a- 1 (d i d j -T k ij (g)d k )a 

- {T k k i(G) + a~%a) {T m mj (G) + a~ l d 3 a) (90) 

+2T k mk (g)v m i:j (g) - v k mj (g)T? k (g) 

~^~9 9 \_Skrsfmij gkm{iQ 'j)rs ^?fcrs^?(ij')m] ■ 

Together with fl28|) , (j83|) and (]89f) constitute a symmetric hyperbolic system 
for the evolution of gij, Kij and f k ij- 

Note that once f k ij or, equivalently, g k ij are introduced as variables, the 
relation (^) becomes an initial condition and does not a priori hold for all 
time. Equation (|83|) can be related to metric compatibility by putting it in 
the form 

^9k(iR 3 )o = B g kij + d k (2NK ij ) (91) 
-Ag k{i NV m {K 3)m - g j)m K) . 

Here, one sees that if the momentum constraint is satisfied for all time, then 

gkij = 2F(ij) k = d k 9ij (92) 

and the connection is metric compatible. If the momentum constraint is vio- 
lated, metric compatibility is sacrificed. This shows the price paid to achieve 
a symmetric hyperbolic system that is close to the canonical equations: the 
momentum constraints become dynamical, and metric compatibility is lost 
if the latter are violated. 



8 Conformal "Thin Sandwich" Data for the 
Initial Value Problem 

The standard approach to the initial value problem is the "conformal method," 
the fundamental rudiments of which were introduced by Lichnerowicz ||. 
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The essentially complete form was developed by two of us (YCB and JWY), 
34 , ||. Basic theorems were obtained by us and by O'Murchadha |45|, 



sec 



5[ [T(| , and Isenberg and Moncrief j|7] . The older approach concentrates (in 
the vacuum case to which we restrict ourselves) on the construction of the 
spatial metric gij = ip A, jij and the traceless part = ifj~ 10 X^ of the extrinsic 
curvature, where = Aij + ^g^K. Here, 7^ is a proper riemannian metric 



given freely, \ %3 is constructed by a tensor decomposition method |34], |48| , 
and K is given freely (not conformally transformed). Note that one may as 
well assume det(7^) = 1 because only the conformal equivalence class of the 
metric matters: the entire method is "conformally covariant." 

N.B.: In Sect. |8], only spatial metrics will be used. Therefore, all overbars 
are dropped in this final section. 

Here, we discuss a new interpretation of the four Einstein vacuum initial- 
value constraints. (The presence of matter would add nothing new to the 
analysis.) Partly in the spirit of a "thin sandwich" viewpoint, this approach 
is based on prescribing the conformal metric [1] on each of two nearby space- 
like hypersurfaces ( "time slices" t = t' and t = t' + St) that make a "thin 
sandwich" (TS). Essential use is made of the understanding of the slicing 
function in general relativity. The new formulation could prove useful both 
conceptually and in practice, as a way to construct initial data in which one 
has a hold on the input data different from that in the currently accepted 
approach. The new approach allows us to derive from its dynamical and 
metrical foundations the important scaling law A % i = ip~ 10 \ l: > for the trace- 
less part of the extrinsic curvature. This rule is simply postulated in the 
one-hyper surf ace approach. 

The constraint equations on S are, in vacuum, 

Vj(K ij -Kg ij ) = , (93) 

R(g) - KijK ij + K 2 = (94) 

where R(g) is the spatial scalar curvature of g^, Vj is the Levi-Civita con- 
nection of g^] and K is the trace of K^, also called the "mean curvature" of 
the slice. 

The time derivative of the spatial metric g^ is related to Kij, N, and the 
shift vector (3 l by 

d t9ij = -2NK l3 + (V l /3 J + V,Pi) , (95) 
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where f3j = gji(3 l . The fixed spatial coordinates x of a point on the "sec- 
ond" hypersurface, as evaluated on the "first" hypersurface, are displaced by 
f3 l (x)5t with respect to those on the first hypersurface, with an orthogonal 
link from the first to the second surface as a fiducial reference: A — ^ * 
where * is the physical spacetime inner product of the indicated natural basis 
four- vectors. The essentially arbitrary direction of J| is why N(x) and (3 l {x) 
appear in the TS formulation. In contrast, the tensor is always deter- 
mined by the behavior of the unit normal on one slice and therefore does not 
possess the kinematical freedom, i.e., the gauge variance, of Therefore, 
N and (3 % do not appear in the one- hypersurface I VP for (T,,g,K). 

Turning now to the conformal metrics in the IVP, we recall that two 
metrics ^ and 7^ are conformally equivalent if and only if there is a scalar 
ip > such that = ip Ar )ij. The conformally invariant representative of the 
entire conformal equivalence class, in three dimensions, is the weight (—2/3) 
unit-determinant "conformal metric" g^j = g~ l ^gij = 7 _1/,3 7ij with g = 
det(<7y) and 7 = det(7y). Note particularly that for any small perturbation, 
g l ^5gij = 0. We will use the important relation 

g ij d t gij = l ij d tgij = gVdtfc = . (96) 

In the following, rather than use the mathematical apparatus associated 
with conformally weighted objects such as g~ij, we find it simpler to use ordi- 
nary scalars and tensors to the same effect. Thus, let the role of g~ij on the 
first surface be played by a given metric 7^ such that the physical metric 
that satisfies the constraints is gij = ip A, jij for some scalar ip > 0. (This 
corresponds to "dressing" the initial unimodular conformal metric g^ with 
the correct determinant factor g 1 ^ 3 = ^ 4 7 1 / 3 . This process does not alter the 
conformal equivalence class of the metric.) The role of the conformal metric 
on the second surface is played by the metric 7^ = 7^ + Uijdt, where, in 
keeping with (|96|), the velocity tensor Uij = d t jij is chosen such that 

. = l ij d tJij = . (97) 

Then, to first order in 8t, 7^ and 7^ have equal determinants, as desired; 
but 7y and 7^ are not in the same conformal equivalence class in general. 

We now examine the relation between the covariant derivative operators 
Di of jij and Vj of g^. The relation is determined by 

r l jk (g) = P ifc ( 7 ) + 2V" 1 {25\ j d k) i, - 7 a 7;*AV>) , (98) 



23 



from which follows the scalar curvature relation first used in an initial- value 
problem by Lichnerowicz [[5], 

R(g) = ^(7) - 8<T 5 A 7 t/> , (99) 

where A 7 = ^D^Diip is the scalar Laplacian associated with 7^. 
Next, we solve (|)5|) for its traceless part 

9t9ij - \gn9 kl dt9ki = Vij = -2NA id + (L g (3) l3 (100) 
with Aij = Kij — (l/3)Kgij and 

(Lgf3)ij = + V,-A - {2/3) 9ij V% . (101) 

Expression (|101|) vanishes, for non- vanishing j3 l , if and only if admits 
a conformal Killing vector (3 l = k l . Clearly, k l would also be a conformal 
Killing vector of 7^, or of any metric conformally equivalent to gij, with no 
scaling of k l . This teaches us that j3 % does not scale. That (3 % does not scale 
also follows because, as generator of a spatial diffeomorphism, it is not a 
dynamical variable. We take the latter "rule" as a matter of principle. 

It is clear in ( |100| ) that the left hand side Uij satisfies ity = ip A Vij because 
the terms in ip cancel out. Furthermore, a straightforward calculation shows 
that 

where V' -4 /?-,- = 7ijV3 l . Next, we note that the lapse function iV has essential 
non-trivial conformal behavior. This is a new element in the IVP analysis. 
The slicing function a(t, x) > can replace the lapse function N, 

N = g 1/2 a . (103) 

(The treatment here extends the one in E{| in a simple but interesting way.) 
We have concluded that a is not a dynamical variable and therefore does not 
scale. Furthermore, without loss of generality, we can set det(jij) = 1 and 
thus g 1 / 2 = ip 6 . Then 

N = tfj 6 a . (104) 

Finally, we fix K and require that it does not scale, as in the standard 
treatment of the IVP [pD| . This step is absolutely essential for geometric 
consistency, as we shall see. 
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Next we solve 



find 



2a 



The momentum constraint becomes 
D 



for A % \ using the scaling rules established above, and 
1 







' 1 ij 




= D 3 


— u 




_2a 



while the Hamiltonian constraint becomes 

8A^ - R{j)$ + ( llkljl )A^A k ^- 7 - (2/3)^V 



. 



(105) 



(106) 



(107) 



The unknowns (ip, (3) obey equations of the same form as do the conformal 
scalar potential and the vector potential W l in the standard analysis |34|. M . 
but no tensor splittings are required. Further, (|106|) and Ql07Q are coupled 
in only one direction when K = const. 

Now we note two interesting consequences of this approach. First we see 
that from iV = ip 6 a, we have identically 



N = g 1/2 a 



(108) 



as a consequence of the method. Therefore, time slices t and t + 5t have a 
relation that is manifestly "harmonic:" 



d N + N 2 K = Nd \oga 



(109) 



a result that is fully consistent with our previous discussions and requiring 
that K be a fixed, non-scaling, variable. 

Finally, we can establish the final relationships between the full rieman- 
nian metrics (t) and g'^ = gij (t + 5t) on the two manifestly harmonically 
related slices t and t + St. As in (|95|), 



dt9ij = dtiip 4 ^) = g ik gji 

Working out ( |110p gives 

9t9%j 



-2N{A kl + -g kl K) + {y k l3 l + V(3 k ) 
3 



t/> 4 [uij + jijd^iplogi/;)) 
Vij + gijd t (ij>]ogi()) , 



;no) 



iii^ 
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where 

WlogV) = 2 -{D k f3 k + Qf3 k d k \ogi,-aK^) 

= d^h) 1 ' 2 = d t {gY" = 2 - {V k (3 k - NK) . (112) 

Hence, d t ip and d t Qij are fully determined and we note that the no-scaling 
rules for f3 k and K were essential. 
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